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Abstract 

The properties of a wide variety of growing models, generically called X/RD, involving the de- 
position of particles according to competitive processes, such that a particle is attached to the 
aggregate with probability p following the mechanisms of a generic model X that provides the cor- 
relations and at random (Random Deposition (RD)) with probability (1 — p), are studied by means 
of numerical simulations and analytic developments. The study comprises the following X mod- 
els: Ballistic Deposition, Random Deposition with Surface Relaxation, Das Sarma-Tamboronea, 
Kim-Kosterlitz, Lai-Das Sarma, Wolf- Villain, Large Curvature, and three additional models that 
arc variants of the Ballistic Deposition model. 

It is shown that after a growing regime, the interface width becomes saturated at a crossover time 
{tx2) that, by fixing the sample size, scales with p according to tx2{p) oc {p > 0), where y 

is an exponent. Also, the interface width at saturation (Wgat) scales as Wsat{p) p~^, {p> 0), 
where S is another exponent. 

It is proved that, in any dimension, the exponents 5 and y obey the following relationship: 
<^ = vPrd, where Pro = 1/2 is the growing exponent for RD. Furthermore, both exponents 
exhibit universality in the p — ^ limit. 

By mapping the behaviour of the average height difference of two neighbouring sites in discrete 
models of type X/ RD and two kinds of random walks, wc have determined the exact value of the 
exponent 5. When the height difference between two neighbouring sites corresponds to a random 
walk that after walking (n) steps returns to a distance from its initial position that is proportional 
to the maximum distance reached (random walk of Type A), one has 5 = 1/2. On the other hand, 
when the height difference between two neighbouring sites corresponds to a random walk that after 
(n) steps moves (/) steps towards the initial position (random walk of Type B), one has 5 = 1. 

Finally, by linking four well-established universality classes (namely Edwards- Wilkinson, Kardar- 
Parisi-Zhang, Linear-MBE and Non-linear-MBE) with the properties of Type A and B of random 
walks, eight different stochastic equations for all the competitive models studied are derived. 

PACS Numbers : 68.35.Ct; 05.40.-a; 02.50.-r; 81.15.Aa 

Key words : Irreversible growth processes; Dynamic Scaling; Interfaces and surfaces. 
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I. INTRODUCTION 



The study and understanding of the properties of growing interfaces have attracted great 
interest. In fact, interfaces are ubiquitous in Nature and their study has opened a promising 
field of multidisciphnary research . Interfaces naturally emerge in a wide variety 

of systems such as film growth by vapour deposition, chemical deposition or molecular 

m " 



beam epitaxy 



propagation of fire fronts 6], diffusion fronts |7 [abacterial growth 



solidification j9l], propagation of reaction fronts in catalysed reactions [10|, electrodeposition/ 



dissolution experiments sedimentation etc. 

Models of growing interfaces may be defined and studied either by using discrete lattices 
or by means of continuous equations. Discrete models are defined by a set of rules that 
provide a detailed microscopic description of the evolution of the growing aggregate. In 
these models the interface is described by a discrete set h{i, t) that represents the height 
of site i at time t. The interface has L'^ sites, where L is the linear size and d is the 
dimensionality of the substrate. The interface of the aggregate is characterised through the 
scaling behaviour of the interface width W{L,t) = \/l/L^y^,f_^\h(i,t) — {h{t))]'^. For this 



purpose, the Family- Vicsek phenomenological scaling approach 



121 has proved to be 



very successful for the description of the dynamic evolution of growing interfaces. In fact, it 



may be expected that W{L,t) would show the spatio-temporal scaling behaviour given by 
Wsat oc L° for tc and W{t) oc for t <^ tc, where tc is the crossover time 



may be < 
HQ: 



between these two regimes. The scaling exponents a, (3 and Z = a/ (3 are called roughness, 
growth and dynamic exponents, respectively. Also, different models can be grouped into 
universality classes when they share the same scaling exponents. 

In contrast to the microscopic details of the growing mechanisms of the interface, contin- 
uous equations focus on the macroscopic aspects of the roughness. Essentially, the aim is to 
follow the evolution of the coarse-grained height function /;,(x, t) by using a well-established 
phenomenological approach that takes all the relevant processes that survive at a coarse- 
grained level into account. This procedure normally leads to stochastic non-linear partial 
differential equations that, in general, may be written as follows [ll Il4 . Il5l. Ilfij] 



g/i(x, t) 
dt 



G,{h{^,t)} + F + r]{^,t), 



(1) 



where the index j symbolically denotes different processes, Gj{h{:K,t)} is a local functional 
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that contains the various surface relaxation phenomena and only depends on the spatial 
derivatives of h{'x, t) since the growth process is assumed to be determined by the local 
properties of the surface only. Also, F denotes the mean deposition rate and 77(x, t) is the 
deposition noise that determines the fluctuations of the incoming flux around its mean value 
F. It is usually assumed that the noise is spatially and temporally uncorrelated. 

In order to establish the correspondence between a continuous growth equation and a 
discrete model one can apply at least three different methods: (i) to numerically simulate 
the model and compare the obtained scaling exponents with those of the corresponding 
continuous equation, (ii) to develop a set of plausibility arguments using physical principles 
and (iii) to derive the continuous equation analytically starting from a given discrete model. 

There are few papers in the direction of the last method. For example, a systematic ap- 



proach proposed by Vvedensky et al 



17[ , where the continuous equations can be constructed 



directly from the growth rules of some discrete models, based on the master equation descrip- 
tion, has been requires a regularisation 
step, in which non-analytic quantities are expanded and replaced by analytic approaches, 
e.g., the step function is approximated by a shifted hyperbolic tangent function expanded in 
a Taylor series. As pointed out by Pfedota and Kotrla the choice of the regularisation 
scheme for the step function is ambiguous. Thus, the coefficients entering in the derived 
continuum stochastic equation cannot be determined uniquely. Another method has shown 
the connection between the ballistic deposition discrete model and the Kardar-Parisi- Zhang 
(KPZ) equation in d = (1 + 1) dimensions. However, this method is not successful in 
d = (2 + 1) dimensions |2l| . 

It is worth mentioning that most of the already mentioned progress in the understanding 
of the properties of interfaces has been achieved when the growth of the aggregate is due 
to one kind of particle only. In contrast, less attention has been drawn to the study of 
the dynamics of competitive growing processes. It is well known that these competitive 
processes are significant to the growth of real materials in at least two different ways: a) 
when the growing process involves two or more kinds of particles, and b) when deposition 
of a single kind of particle is considered, but such type of particle may undergo different 
growing mechanisms. 

One exanaple of case a) arises from the deposition of alloys or systems with impurities, 
see e.g. 

0, Q, Q, Q, 13, Q and references therein. In this case, there may be different 
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interactions between different kinds of particles causing tlie gro wing meclianisms to cliange 
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27|. Based on these ideas, Cerdeira et al j2z 



25| liave stud- 



ied various models for binary systems involving competitive random-like and ballistic-like 
deposition. Recently, the scaling behaviour of a two-component surface-growth model has 
been studied by Kotrla et al 1271]. This study addresses the relationship between kinetic 
roughening and phase ordering in a (1 -|- 1) -dimensional single-step solid-on-solid model with 
Ising-like interactions between two components. 

On the other hand, considering the deposition of one kind of particle (case b), Pellegrini 

nn 

et al [2^ |2^ have studied a ballistic model of surface growth that considers "sticky" and 
"sliding" particles. The model interpolates between a standard ballistic model when only 
sticky particles are deposited (with probability P = 1) and a completely restructured ballistic 
model for P = when only unrestricted slidin g pa rticles are allowed to become attached to 
the sample. Using this model Pellegrini et al |28l. l29 1 have given evidence of a roughening 
transition in dimensions d = 3 and d = 4, while such kind of transition is no longer observed 
md = 2. 

In a related context of competitive growing processes, we have also studied two competi- 
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In the first discrete 



tive growth models in (1 + 1), (2 + 1), and (3 + 1) dimensions |30 . 
growth model, namely the RDSR/RD model, the same types of particles are aggregated ac- 
cording to the rules of random deposition with surface relaxation (RDSR) with probability 
p and according to the rules of random deposition (RD) with probability (1 — p)|2l|- In the 
second discrete growth model, namely the BD/RD model, particles are aggregated accord- 
ing to the rules of ballistic deposition (BD) with pro bability p and according to the rules of 
random deposition (RD) with probability (1 — p)|31|- 

For both the RDSR/RD and the BD/RD models the saturation process of the interface 
width depends sensitively on p: saturation takes place at longer times for smaller values of 
p, while the final width of the interface is smaller for larger p-values. Furthermore, in both 
models, three different regimes and two corresponding crossovers can easily be observed. 
For short times, say t < t^i, the random growth of the interface is observed (i.e. the 
RD process dominates). At this stage, correlations have not been developed yet and W{t) 
oc (t < txi^PRD = 1/2) holds. During an intermediate time regime, say t^i < t < 1^2, 
correlations develop since the RDSR (BD) process now dominates leading to W{t) oc t^^RDSR 
{W{t) oc t^^°). At a later stage, for t > tx2, correlations can no longer grow due to the 
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geometrical constraint of the lattice size and saturation is observed. The saturation value of 



the interface width Wsat{L,p) and the characteristic crossover time tx2 behave as jSO, 
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and 



(p>0), 



(2) 



{p > 0), 



(3) 



respectively. Here, 6 and y are exponents and X = RDSR or BD, depending on the model. 
On the other hand, one has that the crossover time t^i also scales with p as does (see 
equation 

In these previous studies we have shown that the exponents y and 6 are independent of 
the dimensionalit_i^For the RDSR/RD (BD/RD) model we have found that 5 ^ 1(^ 1/2) 

|. Based on these numerical estimates we have conjectured the 



and y 



1 ^ 
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following exact values 6 = 1 and y = 2 for RDSR/RD, and 6 = 1/2 and ?/ = 1 for BD/RD. 
Very recently, this early conjecture has proved to be correct by using an exact analysis 
Furthermore, these values allowed us to formulate another conjecture by statintg that 
6 = y/2 for both models 

Also, we have shown that the stochastic representation of the RDSR/RD model is given 
by 0,0,0 



dh{-K,t) 
dt 



F + Uop'^V^h{x,t) + T]{x,t), 



(4) 



while for the BD/RD model one has 



dh{x, t) u,pV'h{x, t) + ^ [V/i(x, t)]' + r/(x, t), 



(5) 



dt - ■ - ■ V--' - ■ 2 
where plays the role of an effective surface tension, A represents the lateral growth, and 
r7(x, t) is the uncorrelated white-noise term. 

Within this context, the aim of this work is to perform a systematic study of a 
wide variety of competitive growth models of the type X/RD, with X = DST (Das 
Sarma-Tamboronea) 0], KK (Kim-Kosterlitz) 3^, LDS (Lai-Das Sarma) 3^, WV 
(Wolf- Villain) 37|, LC (Large Curvature) 0, RDSR 0, BD |l|, BDl, BD2 and BD3. 
The last three ballistic deposition-motivated models are variants of the BD model that 



will be described in detail below. The study is focused on the behaviour of the exponents 
y and 5 as well as on the derivation of the stochastic equations of competitive models. 
For this purpose the manuscript is organized as follows: in section II and III the relevant 
properties of y and 5 are addressed and the exact relationship between both exponents 
is derived, respectively. Subsequently, in section IV numerical results obtained by means 
of Monte Carlo computer simulations covering all competitive models listed above are 
presented. After that, in section V, the exact values for the exponents y and 5 are 
obtained for all the studied models. Section VI is devoted to the derivation of the stochas- 
tic equations for different competitive models while our conclusions are stated in section VII. 



II. THE RELATIONSHIP BETWEEN 5 AND y 



In previous papers 
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32] we have shown that for the RDSR/RD and the BD/RD 



models, the values of the exponents 5 and y are independent of the dimensionality of the 
substrate. Furthermore, our study leads us to conjecture a simple relationship between them, 
namely 5 = y/2. Furthermore, we have also proposed and numerically tested the following 
)henomenological dynamic scaling Ansatz for both the RDRS /RD and the BD /RD models 



3G 
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W{t,L,p) oc L^^'p-^Fi 



t 



p> 



,t > txl 



oo 



(6) 



where X = RDSR or BD depending on the model and F is a suitable scaling function. Now, 
if we restrict the previous Anzats by considering variations of p only (i.e. fixing L), equation 
(El) can be written as 



W{t,p) (xp-^F*{t/p-y), p>0 ,t>0 ,L^oo, (7) 

where F*{u) is a suitable scaling function such that: (i) F*{u) oc u^rd for u ^ 0, (ii) 
F*{u) oc u^^ for u in the intermediate regimen, and (iii) F*{u) = constant for u ^ 1. It 
is worth mentioning that the Anzats given by equation ((Zj) is valid for the three regimens 
(t > 0) while the previous one, given by equation (jH)), only holds for regimes (ii) and (iii) 
with t > txi- This effect is due to the fact that by fixing the lattice size one also fixes the 
crossover time t^i that depends on L. In order to check the validity of the new Anzats, 
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figure 1 shows log-log plots of W{t, L,p)p^ versus t/p~y as obtained in d = 1 dimensions, 
by using values of p within the range 0.01 < p < 0.64, and for both models. Here, we have 
taken L = 256 and L = 512 for the RDSR/RD and the BD/RD models, respectively. 

Considering the short-time regime t < t^i, or equivalently u — > in equation ((7j), we 
observed that the initial slope is independent of the considered model and it is given by 
Prd = 1/2, since the RD process dominates the early stages of growth. Consequently, 
within the short-time regime, equation ((7j) can be written as 

Wit,p)^p~\t/p-yf'^o^ p>0 ,t<t,, ,L^oo, (8) 

and, according to the results shown in figure 1, this relation has to be independent of p. So, 
this is true only if 



5 = yP, 



(9) 



This result strongly suggests that the factor 1/2 already found in the relationship between 
S and y js^ 31, 3^ is just f3^^. It is also worth mentioning that (3rd is independent of the 
dimensionality of the substrate, and therefore one should expect that equation Q would 
also hold in any dimension. 



III. NUMERICAL EVIDENCE ON NEGLIGIBLE FINITE-SIZE CORRECTIONS 
TO THE VALUES OF THE EXPONENT 6. 

Usually, it is observed that the numerical values of the scaling exponents (a, (3 and 
Z) undergo systematic deviations when the size of the lattices used in the simulations is 
changed. Figure 2 shows log- log plots of Wsat{L,p)/ L°'^° versus p as obtained for lattices 
of different size and taking the BD/RD model. In this figure, we have added the data 
corresponding to L = 10 to the results already shown in paper Using the exact value 
<^BD = ^, straight lines are observed, in agreement with equation (j21), and the best fit of the 
data gives the slope 6 = 0.45 ± 0.01. However, a rather small systematic deviation of the 
data, according to the size of the lattice, is observed: the larger the lattice, the smaller the 
ordinate. This behaviour is due to high-order corrections to scaling that we have neglected 
in equation On the other hand, using the roughness exponent obtained by fitting our 



8 




FIG. 1: Log-log plot of W{t,L,p)/p-^ versus t/p-y where X = RDRS, BD. Results obtained for 
different values of p (0.01 <p< 0.64) and lattices of size L = 256 for X=RDRS and L = 512 
for X=BD. The data corresponding to the BD/RD model have been shifted one decade to the left 
for the sake of clarity. The dotted line has slope f3jiD = 1/2 and has been drawn for the sake of 
comparison. More details in the text. 




FIG. 2: Log-log plots of WsatiL,p)L~"^^ versus p obtained for lattices of different size, as indicated 
in the figure, and assuming ubd = 1/2- The full line has slope 5 = 0.45 and corresponds to the 
best fit of the data. The insert shows the same scaled plot but obtained assuming a^^j = 0.43 
for L = 512, 256, 128 and a'^^ = 0.46 ± 0.05 for L = 10. Again, the fuh line with slope 5 = 0.45 
corresponds to the best fit of the data. 

data, namely = 0.43 ± 0.05 for L = 512, 256, 128, and = 0.46 ± 0.05 for L = 10, it 
is possible to achieve an excellent data collapsing, as shown in the inset of figure 2. In this 
case, the slope obtained by means of a least-squares fit is also 6 = 0.45 ±0.01. Summing up, 
all results shown in figure 2 point out that the scaling Anzats given by equation (j2)) holds 
for the BD/RD model. Furthermore, we would like to emphasise that the systematic shift 
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10"^ 10"^ 

p 

FIG. 3: Log-log plots of WsatiL,p)/L°'^^^'^ versus p obtained for lattices of different size, as 
indicated in the figure, and assuming annsR = 1/2. The full line has slope 5 = 0.97 and corresponds 
to the best fit of the data. 

of the data observed in figure 2 does not affect the slope of the power law, and consequently 
the exponent 6 is almost independent of the lattice size (up to L = 10 in figure 2). 

Figure 3 shows log-log plots of Wsat{L,p)/ L'^^^^'^ versus p as obtained for the RDSR/RD 
model using lattices of different size. In this figure, we have also added the data correspond- 
ing to L = 10 to the results already shown in paper jscl. Using the exact value urdsr = \ 
for L = 512, 256, 128 and anosR = 0.46 for L = 50, straight lines are observed, in agreement 
with equation 0, and the best fit gives the slope 6 = 0.97 ± 0.01. 

So, these results show that the values of 6 are not appreciably affected by systematic 
deviations even when extremely small lattices are used (up to L = 10 in the examples of 
figures 2 and 3). It is worth mentioning that for both models we have arrived at the same 
conclusion in higher dimensions. More explicitly, we have performed numerical simulations 
by using LxL = 6x6in(i=(2 + l) and LxLxL = 6x6x6in(i=(3 + l) and we 
have found 6 ~ 0.45 and 6 ~ 0.97 independently of the dimensionality, for the BD/RD and 
RDSR/RD models, respectively. (These results are not shown here for the sake of space). 

The lack of appreciable finite-size effects in the values of 6 (up to L = 10) is a rather 
surprising result. It should be noticed that a similar behaviour is also found in the exponent 
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(3rd for the random deposition model, and it is due to the lack of correlations between 
different columns of the aggregate. So, although we have not a convincing explanation for 
the behaviour of 6 in our case, we think that it may be related in a non-trivial way to the 
decorrelation induced by the fraction {1—p) of particles that are deposited according to the 
random deposition rules, which are expected to play a relevant role precisely for p ^ 0. 



IV. SYSTEMATIC EVALUATION OF THE EXPONENT 6 FOR A FAMILY OF 
COMPETITIVE MODELS. 

In this section we take advantage of the almost negligible dependence of the values of the 
exponent 6 on the lattice size, as discussed in the previous section, in order to evaluate it 
by using small lattices for a wide family of competitive models called X/RD. This family 
of models is defined such that particles of the same type are aggregated according to the 
rules of a generic discrete model X with probability p and according to the rules of RD with 
probability (1 — p). 

We expect that for this family of models, the dynamics of the RD process would play the 
same role as in the cases of the RDSR/RD and the BD/RD models. That is, RD causes the 
slowing down of correlations among particles. So, we also expect that equations and (jH} 
would be valid, but the values of the exponents a and Z entering in the scaling relationships 
have to be those of the X model that introduces the correlations among particles. Further- 
more, we expect that the general relationship given by equation © will also hold, allowing 
us to evaluate the exponent y after the determination of 6. 

Figure 4 shows log-log plots of Wsat versus p,md = 1, as obtained for the family of models 

amboronea) js^j, KK (Kim-Kosterlitz) 



of the type X/RD where X = DST (Das Sarma- 

LDS (Lai-Das Sarma) ^3^, WV (Wolf- Villain) 37|, LC (Large Curvature) lasj, RDSR 
BD [1], BDl, BD2, BD3. The last three models are variants of the BD model and the 
difference between them is due to the rules used for the sticking of the particles that are 
shown schematically in figure 5. In fact, for the BDl model, when the height of a selected 
site for the deposition of a particle is lower than that of a neighbouring site, the particle 
becomes stuck at half height between the selected site and the highest neighbouring site. For 
the case of the BD2 model, when the height of the site selected to deposit a particle is lower 
than that of a neighbouring site, the particle becomes stuck to the highest neighbouring 
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FIG. 4: Log- log plot of Wgat versus p for the following competitive models in dimension d = 1: 
DST/RD for L = 6, KK/RD for L = 10, LDS/RD for L = 8, WV/RD for L = 8, LC/RD for 
L = 10, RDSR/RD for L = 10 , BD/RD for L = 10 , BDl for L = 10, BD2 for L = 10, BD3 for 
L = 10. Different plots have been shifted vertically for the sake of clarity. 



site but the actual deposition height is decreased by one lattice unit. Finally, for the BD3 
model, when the selected site for the deposition of a particle is lower than a neighbouring 
site, the particle sticks leaving a single hole in the selected column, so after deposition the 
height of the selected site becomes enlarged by two lattice units. (See also figure 5). 

As shown in figure 4, the behaviour of S exhibits two relevant features. On the one 
hand, the obtained values of 6 are independent of the universality class of the competitive 
model that introduces the correlations to the aggregate. So, models belonging to different 
universality classes may have the same value of S, while models belonging to the same 
universality class may have different values of 5. On the other hand, the exponents 5 
obtained in the p — > limit can only assume two possible values that are very close to 
5 = 1/2 and 5 = 1. Further support to this statement is given in figure 6 that shows plots 
of the effective value of the exponent 6 obtained for the family of X/RD models. So, we 
conjectured that in the p — > limit the exact values for 5 should be 1/2 or 1, depending 
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FIG. 5: Schematic view of the deposition of particles in the BDl, BD2 and BD3 models. More 
details in the text. 

on the model. Furthermore, due to the large number of models considered, it would not 
be surprising that the exponent 5 for all models of the type X/RD may assume one of the 
already found values. 



V. ANALYTICAL CALCULATION OF THE EXACT VALUES OF 6 IN THE p^O 
LIMIT. 

In order to evaluate the exact values of 5 we have used three concepts linking the interface 
width (W) to the height difference between two neighbouring sites {h{i) — h{i + 1)) in a 
competitive growth model. 

First, we take into account that from the statistical point of view, it is easy to show 
that if the interface width stops its growth (i.e. by reaching a saturation value), then the 
height difference between two neighbouring sites must stop its growth too. Secondly, we 
considered that if the growing stage of the interface width can be described by a power law 
of the type W{t) oc t^, then one also has that | h{i,t) — h{i + l,t) |(x t^. Finally, if the 
dependence of the interface width on p is given by another power law, namely W{p) oc p'', 
then the relationship | h{i,p) — h{i + l,p) |oc p^ must also hold. It is worth mentioning that 
the inverse relationships are also valid in all cases. 

So, based on these statements we studied the behaviour of | h{i) — h{i + 1) | in order 
to obtain the relevant properties of the interface width. The simplest example for the 
application of the above-discussed concepts is provided by the RD model. Here, the height 
difference between two neighbouring sites due to the deposition of a single particle can 
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FIG. 6: Log- linear plot of S versus p for the following competitive models in dimension d = 1: 
DST/RD for L = 6, KK/RD for L = 10, LDS/RD for L = 8, WV/RD for L = 8, LC/RD for 
L = 10, RDSR/RD for L = 10, BD/RD for L = 10 , BDl for L = 10, BD2 for L = 10, and BD3 
for L = 10. 

only be increased or reduced by one lattice unit. So, this height difference corresponds 
to the displacement of a typical random walk. This fact allows us to conclude that if 
h{i, t — 0) — h{i + l,t — 0) — (as usual), then t) — h{i + 1, t))'^) oc t and using this 
result we conclude that W{t) oc i^/^, in agreement with the fact that Prd — 1/2. 

Considering a competitive growth model, let us call "RD particles" and "X particles" 
those particles that are deposited according to RD rules and X rules, respectively. Since in 
the p ^ limit the deposition probability of RD particles is larger than that of X particles, 
the deposition process in an X/RD model can be thought as cycles involving the deposition 
of rii RD particles followed by a single X particle (rij is the number or RD particles in the 
cycle i). Also, it is clear that from the statistical point of view, the evolution through this 
kind of cycles is equivalent to the evolution of any hypothetic model such that the same 
number (n) of RD particles followed by an X particle is deposited during each cycle, where 
(n) are the average number of RD particles given by 
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(n) = ^n(l-p)>= P. (10) 

n=0 ^ 

Since the statistical evolution of the height difference {h{i) — h{i + 1)) is related to the 
saturation width of the interface according to 



Wsat{p)^^A{iK^)-hi^ + l)f)), (11) 

our aim is to show that the behaviour of Wsat{p) can also be derived from the statistical 
knowledge of {h{i) — h{i + 1)) as a function of the deposition cycles. 

On the one hand, it is straightforward to show that for models such as BD/RD,BD1/RD 
and BD2/RD, the behaviour of h{i) — h{i + 1) in each cycle corresponds to a random 
walk that after (n) steps returns either (i) to its initial position, (ii) to a point placed 
at a distance that is just half of the maximum distance reached from its initial position, 
or (iii) to a neighbouring site of its initial position. So, these types of competitive models 
{h{i) — h{i + l)) correspond to a random walk that in each cycle walks (n) steps and after that 
returns to a site placed at a certain distance (from the starting point) that is proportional to 
the maximum distance reached from its initial position, a being the proportionality constant. 
Hereafter, this kind of random walk is called "Type A". Now, when the number of cycles 
is large enough the saturation of the interface width is expected to occur and consequently 
one has 

oo 

m) - h{z + m = {n) j2 a' = (12) 

where j is the number of cycles. So, using equation (fTUI) it follows that equation (fT^ is 
equivalent to {{h{i) — h{i + 1))^) oc 1/p {p — * 0), so that equation |TT|) gives 

W{t) oc p-^'"^. (13) 

On the other hand, it is easy to show that for the BD3/RD, TDS/RD, LDS/RD, KK/RD, 
WV/RD, LC/RD and RDSR/RD models, the behaviour of {h{i) - h{i + 1)) in each cycle 
corresponds to a random walk that after (n) steps returns to a site placed at a distance (/) 
from the maximum distance reached during the walk (we call this random walk "Type B"). 
So, after the first cycle one has 
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(((M^)-M^ + l))l)^) = (v^-0^ 

then after the second cycle it follows that 



(14) 



(((M.) - hit + = {\l{n) + {^)-if - /)^ (15) 

and so on. The analytic form of {{h{i) — h{i + 1))^) when the number of cycles is large 
enough can be obtained by following a simple procedure: Since after reaching saturation 
{h{i) — h{i + 1)) stops growing, in the next cycle after saturation the increment of {h{i) — 
h{i + 1)) due to the deposition of (n) RD particles is equal to the effect caused by the 
subsequent X particle. As we have already discussed, the effect of the additional X particle 
is just to move towards its initial point / step. So, if before the deposition of the new particle 
one had {{h{i) — h{i + 1))^) = (where c is a constant), then after that deposition one 
already has {{h{i) — h{i + 1))^) = (c — Z)^. Then, the difference between them is given by 
{{h{{} - h{i + l))2)te/ore) - {{h{{} - h{i + l)f)after) ~ 2l{{h{i) - h{i + 1))), and it has to be 
equal to the effect of the (n) RD particles. So, 

{{h{i)~hit + l)))cx{n), (16) 
and using equation one obtains 

W,M^P''- (17) 

Therefore, we conclude that for all models that can be mapped into random walks of 
Type A and B one has that 5 = 1/2 (see equation (fT!?|) ) and 6 = 1 (see equation (fTTj) ). 
respectively. These results provide an independent, more general, confirmation of the exact 
values of the exponents obtained very recently for the RDSR/RD and BD/RD models js^. 
Of course, by using equation Q one can also obtain the exact value of the exponent y for 
any model that can be mapped into the two types of random walks already considered. 

VI. PHENOMENOLOGICAL STOCHASTIC GROWTH EQUATIONS. 

It is well known that the stochastic growth equations describing the BD/RD and the 
RDSR/RD models are the KPZ (Kardar-Parisi-Zhang, see equation (jS])) and the EW 
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(Edwards- Wilkinson, see equation (jH) ) equations, respectively. So, both models belong 
to the same universality class as that of the model that introduces the correlations among 



particles, namely the X model with X = BD and X = RDSR [30|, |31|, |32, |33[. Also, using 
scaling arguments on p and the values of 6 and y we have found that for the BD/RD model 
the parameter p appears in the linear and non-linear terms of the stochastic equation, taking 
the form u(p) = up and \{p) = Xp^^"^ (see equation (0)), while for the RDSR/RD model the 



31 



32, aj- 



parameter p appears as a factor of the form z/(p) = up"^ (see equation Pj))|30. 

The scaling argument on p implies that, assuming the interface /i(x, p, t) to be self-similar, 
on rescaling the coordinate p according to 



p —>■ p = cp, (18) 



and the height according to 



h^h' = c-^h, (19) 



one should obtain an interface that is statistically indistinguishable from the original one. 
Since the interface roughness depends on time t as well, one should have 

t^t = c~H. (20) 

So, for any competitive model belonging to the X/RD family we can propose a stochastic 
growth equation similar to its X model equation where the parameter p appears in the 
linear and non-linear terms. After that, by using the scaling arguments on p and the exact 
values of 5 and it is possible to obtain the exact dependence of the prefactors on p for 
any stochastic equation. This systematic procedure and the results of the previous section 
allows us to conclude that the behaviour of {h{i) — h{i + 1)) quantitatively determines the 
exact dependence of the stochastic equation on p . 

Summing up, by using this procedure we have found that in the cases of BDl/RD, 
BD2/RD and BD/RD models, the KPZ stochastic growth equation is given by 

= F + upV'hi^, t) + ^[VMx, t)r + r^(x, t). (21) 

On the other hand, for the KK/RD and BD3/RD models, the KPZ stochastic growth 
equation is given by 
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= F + W^Mx, t) + ^[VMx, t)r + ry(x, t). (22) 

Also, for the case of LDS/RD, WV/RD, LC/RD and RDSR/RD models, the EW stochas- 
tic growth equation is given by 

= F + V^/i(x, + ry(x, t), (23) 
while for the case of the DST/RD model, the stochastic growth equation is given by 

dh 

^ = -Kp'W'h{r, t) + X^p'W'{Wh{v, t)r + rj{v, t). (24) 

Finally, in Table I we have summarized the list of all possible stochastic equations re- 
sulting for competitive deposition models belonging to four different universality classes 
when the competitive process can be mapped into the two types of random walks already 
considered. 



Table I: Summary of stochastic equations corresponding to models belonging to four 
different Universality Classes such that the competitive process can be represented by 
random walks of Type A and B. More details in the text. 



Universality Class 


Random Walk "A" 


Random Walk" B" 


Edwards- Wilkinson 


a/i(x, t) /dt = UopV^hix, t) 
+^(x,t) 


dh{x,t)/dt = UoP^y^h{x,t) 
+^(x,t) 


Kardar-Parisi-Zhang 


dh{x, t)/dt = upV^h{y., t) 
+Ap3/2[v/i(x,i)]2 + ,7(x,t) 


dh{yi,t)/dt^up'^V^h{x, t) 
+Ap3[V/i(x,t)]2 + 77(x, t) 


Linear-MBE 


9/i(x, t) /dt = -KpV^h(^, t) 

+r/(x,t) 


dh{K, t) /dt = -Kp'^V^h{^, t) 


Non hnear-MBE 


dh{x,t)/dt = -KpV^h{x,t) 
+Aip3/2v2(VMx,t))2 + 7;(x,t) 


dh{x,t)/dt = -Kp^V^h{:>c,t) 
+Aip3v2(V/i(x, t)y + r;(x, t) 
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IV. CONCLUSIONS. 



We have studied a wide family of competitive growth models (generically called X/RD 
models) where particles of the same type are aggregated according to the rules of a generic 
discrete model X with probability p and according to the rules of random deposition (RD) 
with probability (1 — p). 

First, we have focussed our study on the properties of the exponents related to the 
interface width Wsat oc p~^ and the characteristic crossover time oc We have shown 
that both exponents are not independent and one has that the exact relationship 5 — vPrd 
(Prd = 1/2) holds for the BD/RD and RDSR/RD models. However, we expect that the 
above relationship would hold for any competitive growth model of the type X/RD. 

Also, we have found that the values of the exponent S do not significantly depend on the 
finite size of the sample. This property has allowed us to systematically study competitive 
growth models using lattices of modest size. This study shows that S exhibits universality 
and its values, in the limit p — > 0, are restricted to either 5 — 1/2 or 5 — 1, depending on 
the model considered. Furthermore, by using a correspondence between two neighbouring 
sites in the discrete model {{h{i) — h{i + 1))) and two types of random walks, we have 
determined the exact values of the exponent 5. When the height difference between two 
neighbour sites corresponds to a random walk of Type A that in each cycle walks (n) steps 
and after that returns to a point at a distance proportional to its initial position, one has 
5 — 1/2 and consequently y — 1. On the other hand, when the height difference between 
two neighbouring sites corresponds to a random walk of Type B that after (n) steps moves 
towards its initial position / steps, we have found that 6=1 and y = 2. 

Finally, using the exact vahics for the exponents 5 and as well as a scaling argument 
on p, we have derived the stochastic growth equations for the whole family of competitive 
models studied. So, we conclude that the properties of the height difference at saturation 
{h{i) — h{i + 1)) in the discrete model, which determines the behaviour of Wsat, have 
allowed us to quantitatively determine the exact dependence on p of the cocirse-grained 
stochastic equations. 

It is worth mentioning that the derivation of this type of coarse-grained stochastic 
equations, based on the growing rules of the corresponding microscopic models, is an 
interesting challenge in the field of modern Statistical Physics. So, we expect that the 
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relationships between microscopic parameters and stochastic equations obtained and 
discussed in this work will stimulate studies of this field. 
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